U.S. Department of Energy 

Lawrence 

Livermore 

National 

Laboratory 




Preprint 

UCRL-JC-145031 



A Statistical Hot Spot 
Reactive Flow Model for 
Shock Initiation and 
Detonation of Solid High 
Explosives 



A.L. Nichols, C.M. Tarver 



This article was submitted to 

12 th International Detonation Symposium, San Diego, California, 
August 11-16,2002 



July 1 , 2002 



Approved for public release; further dissemination unlimited 



DISCLAIMER 



This document was prepared as an account of work sponsored by an agency of the United States 
Government. Neither the United States Government nor the University of California nor any of their 
employees, makes any warranty, express or implied, or assumes any legal liability or responsibility for 
the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or 
represents that its use would not infringe privately owned rights. Reference herein to any specific 
commercial product, process, or service by trade name, trademark, manufacturer, or otherwise, does not 
necessarily constitute or imply its endorsement, recommendation, or favoring by the United States 
Government or the University of California. The views and opinions of authors expressed herein do not 
necessarily state or reflect those of the United States Government or the University of California, and 
shall not be used for advertising or product endorsement purposes. 

This is a preprint of a paper intended for publication in a journal or proceedings. Since changes may be 
made before publication, this preprint is made available with the understanding that it will not be 
cited or reproduced without the permission of the author. 



This report has been reproduced directly from the best available copy. 

Available electronically at http:/ / www.doe.gov /bridg e 

Available for a processing fee to U.S. Department of Energy 
and its contractors in paper from 
U.S. Department of Energy 
Office of Scientific and Technical Information 
P.O. Box 62 

Oak Ridge, TN 37831-0062 
Telephone: (865) 576-8401 
Facsimile: (865) 576-5728 
E-mail: reports@adonis.osti.gov 

Available for the sale to the public from 
U.S. Department of Commerce 
National Technical Information Service 
5285 Port Royal Road 
Springfield, VA 22161 
Telephone: (800) 553-6847 
Facsimile: (703) 605-6900 
E-mail: orde r s^nlis-f edw Q r l d tgQY 
Online ordering: http: / /www.ntis.gov/ordering.htm 

OR 

Lawrence Livermore National Laboratory 
Technical Information Department's Digital Library 
http: / / www.llnl.gov/tid/Iibrary.html 



A STATISTICAL HOT SPOT REACTIVE FLOW MODEL FOR SHOCK INITIATION 
AND DETONATION OF SOLID HIGH EXPLOSIVES* 



Albert L. Nichols III and Craig M. Tarver 
Lawrence Livermore National Laboratory 
Livermore, CA 94551 



A statistical hot spot reactive flow model for shock initiation and detonation of solid high 
explosives developed in the ALE3D hydrodynamic computer code is presented. This model is 
intended to evolve into a physically correct description of the physical and chemical 
mechanisms that control the onset of shock initiation via hotspot formation, the growth (or 
failure to grow) of these hotspots into the surrounding explosive particles, the rapid transition to 
detonation, and self-sustaining detonation. Mesoscale modeling of the shock compression and 
temperature dependent chemical decomposition of individual explosive particles are currently 
yielding accurate predictions of hot spot formation and the subsequent growth (or failure) of 
these hotspot reactions in the stirrounding grains. For two- and three-dimensional simulations 
of larger scale explosive charges, a statistical hotspot model that averages over thousands of 
individual hotspot dimensions and temperatures and then allows exothermic chemical reactions 
to grow (or fail to grow) due to thermal conduction is required. This paper outlines a first 
approach to constructing a probabilistic hot spot formulation based on the number density of 
potential hotspot sites. These hotspots can then either ignite or die out if they do not exceed 
certain ignition criteria, which are based on physical properties of the explosive particles. The 
growing hot spots spread at bum velocities given by experimentally determined deflagration 
velocity versus pressure relationships. The mathematics and assumptions involved in 
formulating the model and practical examples of its usefulness are given. 



INTRODUCTION 

Phenomenological reactive flow models for the 
shock initiation and detonation of solid high 
explosives, such as the Ignition and Growth model 1,2 
and the Johnson-Tang-Forest (JTF) model 3 have 
been very successful in reproducing most of the 
main features of these reactive flows. The reaction 
rate expressions in these models depend upon the 
average compressions and pressures attained in the 
reacting explosive mixture rather than the local hot 
spot temperatures, which are known to control the 
reaction rates in the preferentially heated regions of 
the explosive charge. Thus there are some 
situations, such as shock desensitization 4 , that are 
not easily treated by these phenomenological 
models. Mesoscale modeling of the shock 
compression 5 and temperature dependent chemical 
decomposition 6 of individual explosive particles are 
currently yielding accurate predictions of hot spot 
formation mechanisms and the subsequent growth 
(or failure) of these hot spot reactions in the 
surrounding grains. For two- and three-dimensional 
simulations of realistic size explosive charges, a 
statistical hot spot model that averages over 
thousands of individual hot spot dimensions and 



temperatures and then allows the chemical reactions 
to grow (or fail to grow) due to thermal conduction 
is required. Some simple statistical hot spot models 
were developed several years ago in one- 
dimensional hydrodynamic codes 7,8 , but practical 
hot spot models in multidimensional codes are just 
beginning to appear. 

The chemical kinetic decomposition models for 
HMX and TATB used in this statistical model have 
been previously used to determine the critical hot 
spot temperatures and dimensions that will react and 
begin to grow 9 . The growth rates of these hot spots 
due only to thermal conduction were also 
determined 10 , and were shown to be too slow to 
account for measured shock initiation times and run 
distances to detonation. Therefore physical 
mechanisms, such as crack nucleation and growth, 
must be causing greatly increased reactive surface 
area as the pressure and temperature increase during 
shock initiation and the transition to detonation. 

In this paper, the formulation of a statistical hot 
spot creation model in the ALE3D hydrodynamic 
code is presented. During shock compression, a 
distribution of hotspot sizes and temperatures 
based on the initial pore size distribution, explosive 






particle size distribution, and the density of the 
explosive charge is created. These hot spots then 
grow to consume neighboring explosive particles or 
fail to grow as thermal conduction lowers their 
temperature before exothermic chemical reaction 
can be completed. The growth rate of the surviving 
hot spots is then modeled by a statistical hotspot 
growth model normalized to experimental high- 
pressure deflagration rate data. The coalescence of 
these expanding hot spots as the temperature and 
pressure rapidly rise everywhere in the decomposing 
explosive charge is also discussed. The extremely 
rapid transition to detonation observed 
experimentally is shown to be similar to a constant 
volume explosion from a critical compression state 
of the unreacted explosive to a state on the reaction 
product Hugoniot close to the Chapman-Jouguet (C- 
J) state. Examples of calculated shock initiation and 
detonation wave propagation for an HMX-based 
explosive are presented. Examples of shock 
initiation effects, such as desensitization and 
changes in sensitivity due to initial temperature and 
particle size distribution variations, that are difficult 
or impossible to compute with pressure and 
compression dependent reactive flow models are 
also presented. 



STATISTICAL HOTSPOT MODEL 



PROBABILISTIC HOTSPOT FORMULATION 
The first phase in constructing the statistical hot 
spot model is the consideration of the distribution of 
those hot spots. First, consider the probability P r that 
single hotspot of radius R will have reacted at a 
given location in a volume V in the explosive. This 
probability is given by: 



P = 



4teR 3 

3V 



(l) 



If there are Nr of these hotspots randomly 
distributed in space, then the probability that a given 
location has not reacted Pnr is simply the product of 
the independent probabilities. Assuming that the hot 
spots are independently located, Pn,- is defined as: 
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Taking the limit where the volume becomes 
large but the hotspot density remains fixed, Eq. (2) 
becomes: 
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Finally, the combined probability of a region 
not having been reacted Pnr(R) is simply the product 
of the probabilities associated with each hotspot 
radius. Therefore the final expression for P^ is: 



/ 
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(4) 



The probability of not yet reacting is simply the 
mass fraction of the reactant in a reactive flow 
formulation. The probabilistic formulation makes it 
easier to consider a variety of different possibilities. 
For example, a similar line of reasoning can be used 
for two-dimensional hotspots (hotlines) and one- 
dimensional hotspots (hotplanes). The latter would 
be useful for the modeling of shear banding as an 
ignition source, for example. If all of such ignition 
mechanisms could be defined, all that would be 
required for a complete hot spot model is to multiply 
their probability functions together. 



HOTSPOT DENSITY MODEL 

The process derived in the previous section 
defines a mechanism for connecting the probability 
that some quantity of explosive has been reacted to 
the density distribution of hotspots. We define the 
probability density of the hotspots as: 



p(R,t)= p A (R,t)+ p D ( R,t ) (5) 
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the density function, respectively. The number of 
hotspots that are active at time t is A (t), and p B (t) 

is the number of hotspots created at time t. In the 
current model, it is assumed that all hotspots active 
at time t have the same rate of death fi( t), that is: 
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where p s (a,o) is the number density of hotspots that 
ignited at time a, and died at time (0. The dirac-delta 
functions are used to define the size of the hotspot 
with the assumption that the initial hotspot size is 8, 
and that it then bums at a bum rate v out from that 
initial spot. The first term p A (R,t) represents the 
population of hotspots of size R that are still 
growing (active) at time t. The second term pn(R,t) 
represents the population of hotspots of size R that 
have stopped growing (died) by time t. It is 
important to remember that even though a hotspot 
may stop burning, the material that has burned 
within that hotspot must still be counted as reacted. 

Let us now define the following projections of 
the density function: 



p s (a,t)= p(t)\ pXa,(D) (9) 

With these projections, it is possible to 
construct a set of differential equations to couple the 
high order unreacted mass fraction with the much 
simpler active hotspot density. 
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The h term is just the negative of the log of the 
probability defined in Eq. (4). The g and f terms are 
the two- and one-dimensional active projections of 



IGNITION MODEL 

In order to complete the set of equations defined 
in the previous section, we must define the rate at 
which hotspots are created. In order to model the 
explosive process, it is necessary to choose an 
ignition model that can encompass a variety of 
phenomena associated with high explosives. We 
begin by defining the initial density of potential 
hotspots pp°. These potential hotspots can be 
anything from defects in the crystal lattice to voids 
in the region between the explosive grains. For the 
current model, we only limit ourselves in that the 
potential hotspot must transform into a roughly 
spherical hotspot. Most postulated hotspot formation 
mechanisms involving void collapse predict that the 
spherical hotspots form upon full collapse 11 . Other 
potential hotspot formation mechanisms, such as 
shear-banding, would transform into roughly planar 
hotspots and thus are not considered in this 
treatment. The shock process compresses these 



potential hotspots. If they are compressed to a 
sufficiently high temperature, they will start to react 
(a hotspot). If the process is too weak, then the 
potential hotspot will be destroyed without creating 
a hotspot. Without such a process, any sufficient 
compression of the explosive would lead to reaction, 
even that from an isostatic press. However, since 
explosives do have strength, there must be sufficient 
force to overcome that internal void strength before 
any changes to the potential hotspot density can 
occur. The following phenomenological ignition 
model captures these features. 



K = 



' AP*(JP—P 0 ) 
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( 11 ) 



Pp=~Pp K 

p B =p p (K-D)H(K-D) 



equal to the value of R for the shock pressure that 
just begins to ignite the explosive. 

A heuristic argument can be used to determine 
p P ° 9 and £. If one assumes that the initial hotspot 
volume will equal the initial void volume and that 
enough hotspots need to be created so that when 
they bum they will consume the entire explosive in 
the reaction zone time % and initial void density p v : 



Pv « 47Z£ 3 p° 3 

, 0 / (13) 
1 = 47t(£+vtJp°/3 

Another heuristic argument can be applied to 
the ignition pre-factor A and P*. It is reasonable to 
asstime that the rate of collapse of the void regions is 
proportional to the product of the void radius and the 
particle velocity u in the shock wave. For relatively 
low pressures, we can Taylor expand the volume 
change using the adiabatic compressibility, and get: 



Here K is the rate of potential hotspot 
transformation, and D is the constant death rate for 
potential hotspots. Po is the ignition rate threshold 
pressure that represents the internal resistance to 
void collapse. To prevent unrealistically large 
collapse rates during numerical pressure spikes, P is 
defined as the saturation pressure. H is the heavy 
side step function, which is zero for all arguments 
less than zero and one for everything else. We 
originally envisioned a compression rate dependent 
ignition rate, but such a rate can be extremely mesh- 
size dependent. More complex ignition models can 
be formulated as this model evolves. 
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where po is the initial density, and c is the reference 
sound speed. Comparing to the ignition definition in 
the limit of small (P-Po), we find 

Am ska (i5) 



PHYSICAL INTERPRETATION OF HOTSPOT 
MODEL PARAMETERS 

We have defined a total of 8 parameters, not 
counting those associated with the equation of state, 
for this model. They are: Po, P*, A, /z, v, D, p P ° 9 and 
£. Po is clearly related to the yield strength of the 
explosive, and so we will use the yield strength in 
our model. The burn velocity v can be 
experimentally determined by any of the standard 
bum rate measurement techniques, such as strand- 
burner and diamond anvil experiments. The value 
of D should be chosen to match explosive shock 
recovery experiments, so that the value of D is set 



Although it might be tempting to use an ignition 
form defined by the particle velocity, the 
complexities associated with making it work under 
the variety initial conditions would be daunting. For 
example, changes in temperature would necessitate a 
change in the reference density within the root sign. 
Without such a change, hot systems could never 
ignite. We can handle the natural curvature that 
comes out of this formulation by an appropriate 
choice of P*. 




EXAMPLE PROBLEMS 

SHOCK INITIATION 

The phenomenological reactive flow model 
Ignition and Growth 1 has been normalized to a great 
deal of shock initiation data. However, for different 
initial temperatures and particle size distributions, it 
has to be reparameterized. The Statistical Hot Spot 
model has the necessary features to calculate such 
initial condition differences. The first example 
problem is to drive a 2 cm wide by 3 cm long by 0.1 
cm deep block of explosive into a stone wall to 
check the two-dimensional shock initiation 
characteristics of the Statistical Hot Spot model. As 
previously mentioned, an unreacted equation of state 
is needed to describe the states attained during shock 
compression. The Jones-Wilkins-Lee (JWL) 
equation of state is used with typical parameters for 
an HMX-based plastic bonded explosive. 

P = A(1 -oVRiV)e' Rlv — B(l{sflR. 2 V)e 1 ^+00 E/V (16) 

where P is pressure, V is relative volume, E is the 
internal energy, co is the Gruneisen coefficient, and 
A, B, Rj, and R 2 are constants. For a typical HMX- 
based plastic bonded explosive, the initial density is 
1.85 g/cm 3 , Ri = 14.1, R 2 = 1.41, ffi = 0.8938, A = 
9522 Mbar, and B = 0.05944 Mbar. This JWL 
equation fits the measured unreacted Hugoniot data 
at low shock pressures and the von Neumann spike 
data at high pressures 12 . The reaction products are 
described by LEOS tables fit to product Hugoniot 
and isentrope states calculated by the CHEETAH 
chemical equilibrium code 13 . Mixtures of unreacted 
explosive and reaction products are assumed to be in 
pressure equilibrium. 

As discussed in the previous section, the 
ignition and growth of reaction model has eight 
parameters: P 0 , P* 5 A, /t, v, D, p P ° 9 and e. P0 is the 
ignition rate threshold pressure, which is related to 
the effective yield strength at high strain rates and is 
set equal to 0.1 GPa for these calculations. P* is the 
saturation pressure below which all hot spots are 
formed and is set equal to 10 GPa. The ignition 
prefactor A is defined by Eq. (15) and, for HMX, is 
approximately 2000 cm-|is/g. The hot spot death 
rate fi is initially set equal to 5 jis" 1 . The constant 
death rate parameter for potential hot spot D is set 
equal to 11.3. The initial number of potential hot 
spot sites pp is assumed to be 1.4 x 10 10 cm' 3 . The 



initial hot spot diameter e is assumed to be 1.5 x 10" 4 
cm or 1.5 microns. Finally the reaction growth rate 
v is assumed to be a function or pressure as 
measured experimentally in strand burners 14 . This 
pressure versus bum rate function is shown in 
Tablet. 

TABLE 1. PRESSURE VERSUS BURN RATE 



PRESSURE (GPa) 


BURN RATE (cm/lis) 


1.0 x Iff 4 


1.016 x Iff 7 


1.151 x Iff 2 


9.068 x Iff 7 


3.842 x 10' 2 


1.763x1 O' 6 


5.957 x 10' 2 


3.599x1 O' 6 


8.077 x 10' 2 


6.988 x 10- 4 


1.164 xlO' 1 


1.016 x Iff 3 


1.0 


4.441 x 10’ 3 



Using these input values, the explosive was 
driven into the stonewall at various velocities to 
study the buildup of the shock wave produced in the 
explosive to a detonation. In agreement with 
numerous experiments and previous models, the 
buildup toward detonation occurred well behind the 
shock front, the transition to detonation was rapid, 
and the resulting detonation wave traveled at the 
correct velocity, 8.8 mm/jis. Table 2 lists the initial 
velocities and the calculated times to complete 
reaction. These reaction times are reasonable for 
HMX-based plastic bonded explosives. Figure 1 
shows the fraction reacted contours for HMX 
decomposition at 2.33 jis for the 0.3 mm/p,s initial 
velocity case. 



TABLE 2. INITIAL CONDITIONS AND TIMES 
TO COMPLETE REACTION 

VELOCITY fmm/qs^ REACTION TIMES (\is) 



0.18 


2.51 


0.30 


2.27 


0.60 


1.70 


1.2 


1.51 


2.4 


0.29 




0.976 



0.114 




FIGURE 1. Fraction Reacted Contours for HMX 
(Red — no reaction; Blue — Complete reaction) at 
2.33 |xs for an initial velocity of 0.3 mm/ps 

The effects of changing various hot spot 
parameters were investigated using an initial particle 
velocity of 0.18 mm/|is or an initial shock pressure 
of approximately 1 GPa. Reducing the hot spot 
death rate parameter p from 5 to 1 reduced the time 
to complete reaction from 2.51 ps to 1.76 ps, and 
increasing p to 10 caused the fraction reaction to 
remain below 0.45. Lowering the constant death 
rate parameter D to 1.13 from 1 1.3 led to complete 
reaction in 2.3 ps instead of 2.51 ps, and raising D 
to 22.6 forced no reaction to occur in 9.3 ps. 

One of the most important effects on shock 
initiation is the initial particle size effect, which 
controls both the initial number of potential hot spot 
sites and the average size of these sites. Each solid 
explosive exhibits a maximum shock sensitivity at 
certain particle size. If the particles are too large, 
there are fewer sites and these sites are easily 
compressed resulting in relatively low hot spot 
temperatures. 6 If the particles are small, there are 
more sites compressed to high temperatures, but 
many of these hot spots lose their heat by 
conduction before they can react and grow. 6 Table 3 
shows the effects of changing the initial number of 
potential hot spots p P by an order of magnitude in 
each direction. Large particle size formulations 

contain fewer sites and thus are less shock sensitive. 

TABLE 3. EFFECT OF HOT SPOT DENSITY 

Sitesfcm' 3 ) Fraction Reacted Time (lis) 

1.4 xlO 9 0.0523 4.09 

1.4 xlO 10 0.954 1.92 



1.4 x 10 n 

Table 4 shows the effects of initial hot spot 
diameter e on the reaction times. Figure 2 shows the 
fraction reacted for the 0.6-micron diameter case at a 
time of 10.52 ps when the shock wave has traveled 
through the 3 cm long charge. The 0.6 and 0.9 
micron results ignite approximately the correct 
amount of fraction reacted for shock initiation 
studies, and are reasonable hot spot sizes for real 
solid explosive formulations. Since the reaction rate 
falls rapidly with decreasing hot spot size and 
number of hot spot sites, the Statistical Hot Spot 
model is predicting the well-known particle size 
effects on shock initiation of solid explosives. 

TABLE 4. EFFECT OF HOT SPOT SIZE 

Diameter (cm) Fraction Reacted Time (us) 



3.0 x 10' 4 


0.979 


0.29 


1.5 x 10* 4 


0.954 


1.92 


9.0 x 10' 5 


0.1792 


5.86 


6.0 x 10' 5 


0.0322 


10.52 



The pressure versus reaction rate function for hot 
spot growth shown in Table 1 is similar to 
deflagration rates observed in strand burners at 
pressures below 1 GPa. For the higher pressures 
reached during shock initiation and detonation, the 
growth rate function was fit to recent diamond anvil 
cell (DAC) bum rate data on pure HMX 15 by 
Reaugh 6 These rates are shown in Table 5 and are 
much faster than those in Table 1 . 




FIGURE 2. HMX Fraction Reacted Contours at 
10.52 ps for 0.6 pm Diameter Hot Spots impacted 
at an initial velocity of 0.18 mm/ps 



TABLE 5. PPESSURE VERSUS DAC BURN 
RATES 

PRESSURE (GPa) Bum Rate (cm/ps) 

1.0 xlO- 4 2.35 xlO’ 7 

1.0 xlO' 1 5.0 x 10‘ 5 

3.0 7.0x1 O' 4 

50 9.0 xlO* 2 

200 0.9852 

With this higher bum rate function, less hot spot 
ignition is required for realistic reaction buildup 
times. Table 6 shows the effects of hot spot size in 
conjunction with the DAC based growth rates. 
Figure 3 shows the fraction reacted contours at 2.97 
ms for the 1 .5 micron hot spot diameter case. 

TABLE 6. HOT SPOT SIZE EFFECTS WITH 
DAC BURN RATES 



Diameter(pm) 


Fraction Reacted 


Time(fis) 


0.9 


1.0 


0.84 


0.6 


1.0 


1.44 


0.3 


1.0 


2.36 


0.15 


1.0 


2.97 


These examples 


of shock initiation 


calculations 



show some of the basic capabilities of the Statistical 
Hot Spot model. The current parameters are of the 
correct order of magnitude and will be fine tuned for 
specific explosive formulations. 




FIGURE 3. HMX Fraction Reacted Contours at 
2.97 |ls for 1.5 |lm diameter hot spots impacted at 
0.18 mm/|ls 



COLLIDING AND DETONATION WAVES 

The other main area where phenomenological 
reactive flow models sometimes have problems is 
when multiple shock or detonation waves collide, 
yielding regions of very high pressure and reaction 
rates. To test the ability of the Statistical Hot Spot 
model to handle such problems, the explosive block 
described above was initiated on the top and bottom 
by imparting initial pressures to start the reaction 
sequence. Figure 4 shows the collision of the curved 
detonation waves with unreacted HMX at the edges 
still being consumed. Since the reaction rates at 
very high pressures (essentially twice the detonation 
pressure in this example) can be limited in the 
Statistical Hot Spot model, the calculation times do 
not increase dramatically as they can using models 
with pressure dependent reaction rates. The ALE3D 
code can overcome the deformed zone problems that 
purely Lagrange codes often cannot under these 
conditions. This capability will be extremely 
important when the complex three-dimensional 
Mach stem structure of real detonation waves is 
simulated in 3D. 

One of the most important and complex 
modeling problems in solid explosives is the 
phenomena of shock desensitization or dead 
pressing. The most famous experimental example 
is that of Campbell and Travis 4 , in which a PBX 
9404 is detonated at one end and a weak shock wave 
from an underwater detonation impacts the same 
PBX 9404 charges at the other end. If this shock 
wave is too weak to collapse all the potential hot 




FIGURE 4. Collision of two HMX detonation 
waves produced by 5 GPa input pressures 




spot sites, then the detonation propagates through 
the precompressed material. If the shock wave is 
strong enough to start hot spot ignition and growth, 
then the detonation wave continues through the 
partly reacting explosive. However, if the shock 
wave pressure is within a certain range (0.7 to 2.4 
GPa for PBX 9404) which collapses all the potential 
hot spots without causing too much reaction growth, 
then the detonation wave will fail to propagate when 
it reaches a depth at there are only dead or burned 
out hot spots and unreacted particles. Composition 
B, 4 the TATB-based explosive LX-17 16 , and other 
explosives have also been shown to exhibit this 
phenomena. Phenomenological reactive flow 
models can be parameterized to not allow reaction in 
certain compression regimes 16 or to use reaction rate 
limitations , but they can not predict the times 
required for the death of a detonation wave. 
Figure 5 shows the collision of a detonation wave 
with a region of explosive compacted by a shock 
pressure of 1 GPa and the resulting failure of the 
detonation wave to cause further reaction in the pre- 
compressed explosive whose hot spot sites have 
already died or reacted. Therefore the Statistical 
Hot Spot model is definitely capable of 




Figure 5. Shock Desensitization Example. 
Collision between a shock and a detonation at 
0.9 |ls. Figure on left is the extent of reaction, 
pressure is on the right. Note the pressure 
wave continues after the collision but the 
extent of reaction does not. 



quantitatively modeling shock desensitization and 
other pre-conditioning effects on shock initiation 
and detonation wave propagation in three 
dimensions. 

CONCLUSIONS 

In this paper, a statistical hot spot model was 
derived that can be applied to a number of situations 
that are normally not possible to model with a more 
standard reactive flow treatment. Examples were 
shown for shock initiation and colliding waves that 
demonstrated the versatility of the model. A great 
deal more research is necessary to study the 
combined effects of the various parameters and to 
make each part of the model as physically realistic 
as possible. 
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